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We discuss Ward identities for strongly interacting fermion systems described by Eliashberg- 
type theories. We show that Ward identities are not in conflict with Migdal theorem and derive 
diagrammatically Ward identity for a charge vertex both in a Fermi liquid, and when a Fermi liquid 
is destroyed at quantum criticality. We argue that Ward identity for the spin vertex cannot be 
obtained within Eliashberg theory. 
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I. INTRODUCTION 

Conservation laws set powerful constraints on the 
forms of spin and charge response functions of interacting 
electron systems. The most useful constraints are asso- 
ciated with the conservation of the the total charge and 
the total spin of fermions, e to t an d Stot- The time inde- 
pendence of the total charge and spin implies that the 
charge and spin susceptibilities, Xcil,^) an d x s (q, fi), 
respectively, should vanish at q = and any finite O. 
This requirement imposes specific relations between self- 
energy and vertex corrections, known as Ward identi- 

ties mi g. 

The subject of this paper are Ward identities for 
strongly interacting fermionic systems in which the self- 
energy is large, but is local - it depends on frequency, but 
not on momentum, i.e. S(fc,w) w E(a;). These theories 
have recently been introduced, both phenomenologically 
and microscopically, in context of quantum criticality in 
high T c and heavy fermion materials Examples in- 
clude marginal Fermi liquid theory [f| , gauge theories [t| , 
and critical theories of CDW 01 an d SDW 0] transi- 
tions with the dynamical exponent z > 1 . In the last 
two cases, the locality of the problem is not imposed by 
the original Hamiltonian, but emerges at criticality due 
to the fact that for z > 1, collective spin or charge bosonic 
excitations are slow modes compared to electrons. The 
locality of the self-energy was first discovered by Eliash- 
berg in the analysis of electron-phonon problem [l2T | , and 
low-energy electron-boson theories with E(/c,w) ks E(ot>) 
bear his name. The frequency-only dependent self-energy 
also emerges in infinite D theories of metal-insulator 
transition 7] and in local theories of heavy- fermion quan- 
tum criticality [jj. However in these cases, the interaction 
is not small compared to a fermionic bandwidth, and our 
analysis will not be applicable. 

The issue of Ward identities for Eliashberg theories is 
somewhat non-trivial as these theories are justified by 
Migdal theorem that states that vertex corrections ST /T 
should generally be small compared to E(w)/w. Mean- 
while, Ward identities imply that Ar and E(k>) scale with 
each other. We show, however, that Migdal theorem 
and Ward identities just probe different limiting forms 
of the full, momentum and frequency dependent ver- 



tex V(q, Q). Migdal theorem states that for VFq >> 0, 
ST/T <C E(ti>)/cj. Ward identity is applicable in the op- 
posite limit of vpq <C f2, and implies that in this limit, 
5T/T - E(w)/w. 

We derive in this paper the Ward identity for the 
charge vertex near QCP. We show that Ward identity 
holds both away from QCP, when the system is in the 
Fermi-liquid regime, and right at QCP, where the Fermi- 
liquid behavior may be destroyed by quantum fluctua- 
tions. We demonstrate that in the latter CcLSG, el deriva- 
tion of Ward identity by explicit summation of diagrams 
is not straightforward and requires some care. 

We also consider the SU (2) Ward identity for the spin 
vertex. We argue that it holds near a CDW instability, 
but cannot be obtained within Eliashberg theory near 
a magnetic QCP. The argument is that the SU (2) Ward 
identity is related to the conservation of the total electron 
spin. Meanwhile, Eliashberg theory near a SDW transi- 
tion is based on the spin-fcrmion model in which a soft 
collective bosonic mode in the spin channel is treated as 
a separate spin degree of freedom which can flip an elec- 
tron spin. We argue that this effective two-fluid model is 
justified only when the momenta of the collective bosonic 
mode is much larger than its frequency: VFq » ^- In 
this limit, Migdal theorem is valid. However, in the op- 
posite limit vpq « Q, relevant to the SU{2) Ward iden- 
tity, the derivation of the two-fluid model near a magnetic 
QCP breaks down, i.e., spin- fermion model is just inap- 
plicable. 

The structure of the paper is the following. We first 
discuss in the SeclTTl how Ward identities are related to 
the conservation laws. In Sec IIIII we demonstrate that 
Ward identities are not in conflict with Migdal theorem. 
In Sec. II VI we derive Ward identity for the charge vertex 
and discuss the difference between a Fermi liquid and a 
non-Fermi liquid. In Sec. |V|we discuss in some detail 
Ward identity for the spin vertex. 

II. WARD IDENTITIES AND THE 
CONSERVATION LAWS 

For E(fc,w) = E(cj), the full fermion-boson vertex T 
depends on incoming and outgoing fermionic frequencies 
lj and Q + lj (Q is a bosonic frequency) , and on a bosonic 
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FIG. 1: The full particle-hole polarization bubble. Solid 
lines - full Green's functions, black triangle - the full vertex 

r(w,n). 



momentum q. Its dependence on fermionic momenta is 
small by Migdal condition and can be neglected. The rel- 
evant Ward identity relates the self-energy with T(q, u>, fl) 
at zero bosonic momentum o: T(q = 0, u>, f2) = T(u>, fl). 
The identity implies that 



£(w + fl) -E(w) 

n 



(i) 



where E(a;) = w + E(u>) We use the normalization in 
which r = 1 for non- interacting fcrmions. The role 
of Q hi enforcing the conservation laws becomes clear 
when one considers a fully renormalized particle-hole po- 
larization bubble n(<7, fl) at q = 0. Both spin and 
charge susceptibilities scale with H(q, fl). In the RPA ap- 
proximation, Xs{q, fl) = 2U(q, fl)/(l — E7II(q', fl)), while 
Xcfefi) = 2n(g,0)/(l + £71% 0)). Hence, if 11(0, fl) 
vanishes, Xs(0, fl) and Xc(0, fl) also vanish. 

The proof that 11(0, fl m ) vanishes for Eliashberg-type 
theories has not been presented in the literature, so we 



discuss the computation in some detail. The full polar- 
ization bubble in Matsubara frequencies is given by the 
diagram in FigQ 



n(o,fi m ) = - 



d pduo„ 



(2n) D + 1 

T(uj m ,il m )G(p,uj m )G(p,uj m + n m ) (2) 

where the fermionic Green's function G(p,io m ) is given 
by 



G(p,u m ) = (iE(u}f, 



(3) 



Linearizing, as usual, the dispersion near the Fermi sur- 
face as e p = vp(p — pp), and replacing the momentum 
integral J d D p/(2ir) D by Nq J de p , we obtain from 



u(o,n m ) = -N J ^r(uj m ,n r . 



w 



de r , 



w (c p - iS(u m ))(e p - i^(uj m + fl m )) 



(4) 



where W is of order of fermionic bandwidth. It is com- 
mon to set W infinite in continuous theories. We will see, 
however, that it is important to keep W finite at inter- 
mediate stages of calculations, and set W — oo only at 
the very end. The need for this procedure follows from 
the fact that at large frequencies H(u) m ) « ui m , T 1, 
and the 2D integral over de p duj m formally diverges. A 
finite W provides the physical regularization of the diver- 
gence as the linearization of the dispersion is only valid 
at vf(p—Pf) < Ep- The frequency integral, on the other 
hand, has to be evaluated in infinite limits. 
Integrating explicitly over e p in Q we obtain 



n(o, n ro ) = iN Q 



du>r, 



2vr T,(oj m + fi m ) - E(w m ) 

r 



log 



(w - it{u m + n m ))(-w - it(w m )) 
(-w - i±{u: m + n m ))(w - it(u m )) 



(5) 



One can easily make sure that at W — > oo, the frequency 
integral in JSJ is the sum of two terms: 11(0, fl TO ) = 
IIi(0, fl m ) + 112(0, fl m ). The first contribution comes 
from large u> rn 3> fl m for which E(w) ~ W. Expand- 
ing in T,(lo + fl) — E(w) in the logarithm, we obtain 



Hi (0, fir 



TF 



do; 



4iv / r^fi™)- 

/O VF 2 + E 2 (w m ) 27T 



(6) 



The integrand is 0{l/W) up to E(ct>) - W, and 
0(VF/E 2 (w)) at larger frequencies. This implies that 
the integral is confined to uj, for which E(w) ~ VF. 
At these frequencies both the self-energy and the ver- 
tex renormalization are irrelevant, i.e., T(u> m ,Q m ) = 1, 



and E(u; m ) = Lu m . Performing the frequency integration 
in iJBJ we then obtain 



ni(o,n r . 



4N 



W duj 



W 2 



2vr 



= N„ 



(7) 



Another contribution comes from small frequencies u> m ~ 
fl m , when uo m and uj m + fl m have different signs, and the 
argument of the logarithm is iir (sgn(uj + fl) —sgnuS). The 
real part of the logarithm is 0(1/ W) at small frequencies 
and can be neglected. Replacing the logarithm in JSJ by 
its argument, we obtain 



n a (o,n m ) = -jv / dw^, 

_n m E(w, 



r(^m? flm) 



fl r 



E(w m ) 



(8) 
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FIG. 2: The lowest-order diagrams for the fermionic self- 
energy and fermion-boson vertex. Wavy line is the bosonic 
propagator xi( w )- 



For definiteness, we set il m > 0. 
Eq. (JTJ, we immediately obtain 



Using Ward identity, 



-N 



duj r . 



= -N 



(9) 



Comparing (JJJ and JJjJ, we see that the low frequency 
contribution, for which Ward identity is crucial, cancels 
the one from high frequencies, and, as a result, 11(0, fi m ) 
vanishes, as it indeed should. 



III. DIAGRAMMATIC DERIVATION OF WARD 
IDENTITY 

A. the relation to Migdal theorem 

A more subtle issue is how to prove Eq. Q dia- 
grammatically. A conventional recipe is to relate the 
equations for the full vertex and the fermionic self- 
energy 0, 0- At a first glance, this would contra- 
dict Migdal theorem that states that in Eliashberg-type 
theories, vertex corrections are much smaller than the 
self-energy. 0|. We will demonstrate, however, that 
Migdal theorem is only valid for vpq 3> while Eq. 

is valid for VFq <C To understand the impor- 
tance of the VFq/Q ratio, note that every Eliashberg- 
type theory can be viewed as the interaction between 
fermions and the "local" spin or charge bosonic propa- 
gator = Jdq d - l xM H El El El The 

lowest order self-energy and vertex corrections due to 
electron-boson interaction are shown in Fig. [21 For the 
self-energy, we have ^3 

E(w m ) oc g 2 N I de k j dQ m G (k, J m + 0, m )xi(^m) 

(10) 

mN sgnuj m (11) 



Evaluating the momentum integral using 

de k 



de k G Q (k,uj m ) 



iu m - e k 



we obtain 



E(w m )ocg / dQ m Xi{Q m ) 



(12) 



We see that the frequency integral is confined to inter- 
nal frequencies which are smaller than the external one. 
When Lo m is small, relevant tt m are also small, and away 
from the critical point X;(0 TO ) can be approximated by 
XI (0). Then 



E(w m ) = Xuj Ti 



(13) 



where A oc g 2 N xi{0). When A > 1, dZ(uj)/duj = A is 
not small. Note that Eq. (|12J) does not change if we 
re-evaluate the self-energy using the full Green's func- 
tion G(k 1 uj) = (iTi(u>) — e/c) -1 for intermediate fermion. 
This follows from the observation that the relation 111|) 
remains valid for full G. 

For vertex correction, when both q and Q are nonzero, 
we have 

ST(q, ui m , Q m ) cx ig 2 J J d D kdiv' m x 

Go{k,uj' m )G (k + q,uj' m + n m )xi(uJ m ~ w' m ) (14) 

This integral is ultraviolet convergent, and hence momen- 
tum integral can be extended to infinity. Integrating over 
momentum first, we obtain 



5T(q,u m , fi m ) cx 



9 2 N 



y/w m + ( VFq y 



We see that the integral over internal frequency uj' m is 
again confined to \it)' m \ < £l m - Assuming that both Q m 
and w m are small, we again can approximate Xi(w m — a/) 
by Xi(0) and obtain 



ST(q,UJ m ,ft m ) cx A 



fir. 



y/W m + (v F qf 



(16) 



We see that in the limit of vanishing q, AT ~ A i.e., 
5T ~ dH(oj)/duj. This is consistent with Ward identity, 
Eq. QJ. However, in the opposite limit vpq 3> fi, <5r is 
reduced by fl/vpq- For on-shcll bosons, il = uq, where u 
is the effective "sound" velocity [2lJ Hence, Ar ~ Xu/vp- 
The ratio u/vf is a small parameter in Eliashberg-type 
theories H Il6l| , and the strong coupling limit of these 
theories implies that A > 1, while Au/wp C 1. In this 
limit, the vertex for the scattering of fermions by on-shell 
bosons is small, in agreement with Migdal theorem. 



B. the selection of the diagrams 

The next issue is how to select series of diagrams for 
5T(q,u> m , fi m ). A simple experimentation shows that for 
Xu/vf <C 1, relevant diagrams for vertex renormalization 
form ladder series in which each fermionic propagator is 
the full one. All non-ladder diagrams are small in Xu/vf- 
To demonstrate this, compare the two second-order di- 
agrams for 8T(q,uj mi n m ) - the ladder diagram and the 
crossed diagram (see Fig|3]a, b). For both diagrams, we 
consider the limit when both q and f2 m are small, but 
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FIG. 3: The second-order diagrams for the boson- fermion 
vertex V(ui,Q): (a) ladder diagram, (b) crossed diagram, (b) 
the diagram with the renormalization of the side vertex. For 
Eliashberg theories, the ladder diagram is much larger than 
the other two. 



vpq <C £l m - The ladder diagram ST L can be straight- 
forwardly calculated in the same way as above and yields 
ST^ ~ A 2 . For the crossed diagram STq\ we have, set- 
ting external u) m = and |k e;ct | = hp, 



ccg 4 d D kd D pdoj' m du^G (k,uj! m )Go()i+ci,uj' m +£l m ) G (p, w^)G (k e2;t +k-p, J m -J^) Xi^+fi-wJl) Xi{^'L) 

(17) 



The integration over k and ui' m proceeds in the same way 
as before and yields a factor of A. The integral comes 
from vanishingly small uj' < f2 m , hence ui' can be ne- 
glected in the rest of the diagram. The remaining integral 
then reduces to 

2 /> 

ST™ oc A-^ / d D pd9du" 
C Xl(0) J 

xGo(p,OG (l-p,-0 x?«) (18) 

where 1 = k ext + k is a vector whose length is 2kp cos 0. 
To avoid large denominators in the integrand in (|18H . 
both p and 1 p must be near the Fermi surface. For 
each given 1, this selects special hot regions on the Fermi 
surface. The Fermi velocities at p and 1 — p in these 
two regions are generically not aligned, hence, in eval- 
uating the momentum integral in l|18|) . one has to inte- 
grate independently over e p and ei_ p . Using J d D p oc 
(Nq/W) J J de p G?ei„ p and integrating each of the two 
Green's functions in (|18|l over its e we obtain 

flf r >«A 2 x?(0) J (/«)) 2 (19) 

where we introduced f(ui) — Xli^) / Xl{^) i = 1- 

The frequency integral converges at u>" ~ upp such 
that the frequency integral gives upf/W ~ ujvp. Sub- 
stituting this into (|19fl . we find that SV§t ~ \[\(u/vf)], 
i.e., it is small compared to ladder diagrams. One can 
verify that the same smallness in Xu/vf emerges when 
one includes the renormalization of the side vertices in 
the ladder series (see Fig|3]c). 

To prove Eq. Q we now need to go beyond estimates 
and explicitly sum up ladder series of diagrams. At this 
stage, we consider charge and spin vertices separately. 




FIG. 4: The series of ladder diagrams for the full fermion- 
boson vertex. The diagrammatic series reduce to the integral 
equation for F. 

IV. THE CHARGE VERTEX 
A. general proof 

We begin with the charge vertex T c (q = 0). The 
bosonic x;(w) may originate either in the charge or in the 
spin channel. In both cases, we define A via S(o; m ) = Xlo 
at small uj. For charge fluctuations, A = g 2 N xi(0), for 
spin fluctuations, A = 3g 2 Noxi(0). Evaluating the low- 
est order perturbative correction to the charge vertex, we 
find that at small fi m and uj mi T c (u m , £l m ) = 1 + A, i.e., 
vertex renormalization contains exactly the same A inde- 
pendently on whether x/( w ) is °f spin or charge origin. 

Further, the ladder series of of vertex correction dia- 
grams reduce to the integral equation for T c (uj rni Q m ) in 
the form 

r c (^ m , n m ) = l + 1\ J J f{u m - lj'j x 

^cKiAj ^20) 

(«£«J - e k )(iE{n m + uj' m ) - e k ) 

We show this graphically in Fig. 0] 
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A similar equation can be obtained for the self-energy. 
Explicitly computing the prefactor in l|15|) and using the 
fact that the self-energy given by (|15[) does not change if 
we replace Go by full G, as long as £ depends only on 
frequency, we obtain 



S(w m + f2 m ) - E(u m ) = fl m + i\ 

E(q» +<)-£«) 

- e fe )(ffi(fi m + uj' m ) - e k ) 



/(u, ro -a4)(21) 



Comparing Eqs. (|21|l and H20|) . we see that they are 
equivalent if T c (uj m , fl m ) = (E(w m + fi ro ) - T,(u m ))/Q m , 
i.e., when S(w m ) and r c (o; m ,f2 m ) are related by Ward 
identity - Eq. 



B. perturbation theory order by order 

It is also instructive to analyze how Eq. Q is recovered 
by summing up vertex correction diagrams explicitly, or- 
der by order, and comparing the result with the self- 
energy. Order by order consideration was used to prove 
Ward identity for electron-phonon systems Q. There, 
series of vertex corrections are geometrical and one can 
straightforwardly sum them up. We will see that in our 
case, the series of vertex corrections are geometrical in 
the Fermi liquid regime, but become non-geometrical at 
the QCP, where A = oo. In the last case, explicit order 
by order summation of perturbation series is not straight- 
forward and requires some efforts. 



1. Fermi liquid 

In a Fermi liquid, f(u>) ~ 1, and = u> m (l + A). 

Then the leading vertex correction calculated with the 
full fcrmionic propagators is A/(l + A). One can easily 



make sure that higher-order vertex corrections form sim- 
ple geometrical series in A/(l + A). Adding them up, we 
immediately obtain 



r c (g = o) = ^ + / A 



1 + A = 



1 + A \1 + X 
S(w + fi) - E(w) 



l+A 



(22) 



2. quantum- critical point 

Suppose now that the system approaches a QCP at 
which the bosonic mode becomes massless. If the spatial 
dimension is below the critical one d cr , the local suscepti- 
bility Xz(0) diverges at QCP what implies that A— > oo 
(d cr = 3 for Ornstein-Zernike form of x(g,0) [Til Hsjp. 
Suppose further that at QCP, the local susceptibility be- 
haves as Xi( w m) ^m 7 ! an d 7 > 0. Examples include 
3D antifcrromagnetic and ferromagnetic QCP [Tsl H3 | 
and marginal Fermi liquid theory where 7 = 0+ 
(i.e., \i (uJrr,,) oc logo;), 2D antifcrromagnetic QCP where 
7 = 1/2 [ll El, ferromagnetic QCP d [H E H3 
and gauge theory 0, where 7 = 1/3, and correspond- 
ing CDW instabilities with the same 7 [Tol IT^| . In all 
density-wave cases, upper critical dimension is D cr = 3. 

For xi(w m ) oc w~ 7 , self-energy acquires a non- Fermi 
liquid form 



Y,(w Tl 



l— "i 7 



0- 



(23) 



where ojq is the normalization factor. Ward identity im- 
plies that the fully renormalized vertex should behave as 
T cx uo^ 1 , i.e., it should diverge at u) = 0. 

Substituting this self-energy into lj2*U|) and evaluating 
vertex corrections iteratively, order by order, we obtain 
for Vl m ,uj m > 0, 



T c (u) m ,Q m ) = 1 + (1 - 7) 



dx 



o (x+^-V S(x) 



l + (l-7) 



where S(x) = (fl m /w ) 7 + x 1 ^ + (1 - x) 1 ^ . Evaluating 
the first few terms, we obtain for, e.g., 7=1/2 and 

Uj'm f^m — 0. 



r c (0,0) = 1 + 0.7854 + 0.7578- 



(25) 



Clearly, the series are not geometrical, and first few terms 
do not give a hint what the full r c (0, 0) actually diverges. 
This implies that the direct, order by order summation 
of vertex correction diagrams is useless at the QCP. 
Despite that order by order summation fails, infinite 



dy 



1 



o \x-y\-y S(y) 



l + (l-7) 



dz 



o \y-z\-y S(z) 



(24) 



r 



series in (12*4^) can indeed be evaluated explicitly. For this 
we introduce 4>(x) via 



T c (uj mi tt m ) = 1 + (1 - 7) 



lo (x + rC^Six) 

and observe that (j)(x) obeys the integral equation 

dy <t>(y) 



(26) 



(x) = 1 + (1 - 7) 



\x-y\y S(y) 



(27) 
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One can verify that the solution of (|27(l is <fi(x) — 
S(x)(u) /QmV ■ Substituting this cj)(x) into iffi?)) and per- 
forming the integration, we obtain 



ujo_ 

\l<rr}, 



n r 



1-7 



(28) 



This is precisely the same result as one would obtain by 
substituting the self-energy Eq. 123fl . into Eq. In 
particular, when £l m tends to zero, and u> m ~ 
full charge vertex diverges as Sl^ 7 . 



V. THE SPIN VERTEX 

Finally, we discuss the spin vertex T s . Here the situa- 
tion is more involved. When xi originates in the charge 
channel, as e.g., near CDW instability , the spin fac- 
tors in the self-energy and the vertex correction diagrams 
are the same (both equal to one), and Ward identity for 
T s readily follows from the summation of the ladder series 
of diagrams. However, if xi originates in the spin chan- 
nel, as near a SDW instability [Tl| , the interplay between 
the self-energy and the vertex corrections depends on the 
spatial anisotropy of xi ■ If the SDW transition is of Ising 
type, i.e., only the z— component of xi is relevant, the 
spin summation yields the same factor 1 both for the self- 
energy and the vertex pH , IbH Il5l ] , and the ladder summa- 
tion recovers Eq. 2] However, if the transition falls into 
Heisenberg universality class, the summation over spin 
components in the first self-energy diagram yields a factor 

8 <5a,8<JB,a — 3, while the summation over spin compo- 
nents in the first vertex correction diagram yields a factor 
(-1) as J2 a ,3 a a,B^I3n^s,a = -°~§ a 0- As a result, the 
analogy between self-energy and vertex correction dia- 
grams is lost for r s . In a Fermi liquid, ladder diagrams 
for r s now form series in — (1/3)A/(1 + A), and the sum- 
mation of the ladder series yields T s = 3(1 + A)/(3 + 4A) 
what is different from 1 + A expected from Eq. 

The fact that Eq. Q for T s is not recovered in the 
summation of the ladder diagrams does not imply that 
the Ward identity based on the conservation of the total 
S z is invalid. Rather, the non-equivalence between self- 
energy and vertex corrections is the consequence of the 
fact that for spin vertex, Eliashberg theory near the SDW 
instability is only applicable for vfQ 3> but not at 

S> VFq where Ward identity must be valid. Indeed, 
the point of departure for the SDW Eliashberg theory is 
the spin-fermion model 0, ^| . This model assumes that 
the low-energy physics is governed by the vector spin-spin 
interaction 



Uint = g^2s(q)S(q) 



(29) 



between a fermion with a spin s(q) = J^k c ! a <7a - 3Ck +<]- B 
and a spin collective mode described by a bosonic field 
S(q). In other words, a static collective mode is formally 



treated as a separate spin degree of freedom. Within 
this description, the total spin of the system is Stot(q) = 
s(q) + S(q). This total spin is indeed a conserved quan- 
tity However, the Ward identity, Eq. Q), involves only 
the electron spin as is evident from its relation to the 
electron particle-hole polarization bubble n(g, £1). If the 
interaction with the collective mode is of Ising type, there 
are no spin-flip processes between s z (q) and S z (q), and 
s z (q = 0) and S z (q = 0) are two independent conserved 
quantities. In this situation, Eq. is satisfied, as we 
have demonstrated above. However, for Heisenberg-type 
exchange with the collective mode, only Sj ot (0) is con- 
served, but not s z (0) as the electron spin can flip under 
the action of S + (q)s~(q) and transfer its z component to 
S(q). As a result, within the spin-fermion model, there 
is no conservation of the electron spin. 

In reality, S is indeed an auxiliary field. Collective 
modes are made out of electrons, so there is no indepen- 
dent spin degrees of freedom other than electron spins. 
The physics behind the effective spin-fermion model with 
two spin degrees of freedom is the separation of scales. 
In perturbation series, the static part of the spin suscep- 
tibility comes predominantly from electrons with ener- 
gies comparable to the fermionic bandwidth W ~ vfPf, 
while the dynamics of the spin response function is domi- 
nated by the Landau damping and comes from fermions, 
with energies comparable to the frequency of the collec- 
tive mode ft. Eq. I|29|) is valid if relevant f2 are much 
smaller than W. The largest bosonic momenta are of or- 
der pf, hence the largest f2 for typical, mass shell bosons 
is of order upf ■ The condition Q << W then implies 
u << vf, which is precisely the applicability condition 
for Eliashberg theory. This condition, however, also im- 
plies that at small q, typical f2 ~ uq vpq, i.e., the limit 
of q = and finite f2, relevant for the conservation laws, 
cannot be reached. 



VI. CONCLUSIONS 

To summarize, in this paper we considered Ward iden- 
tities for charge and spin vertices in strongly coupled 
Eliashberg- type theories in which T,(k,ui) w £(oj). We 
explicitly demonstrated how Ward identities impose con- 
servation laws and argued that Ward identities are not 
in conflict with Migdal theorem. For the charge vertex, 
we demonstrated that Ward identity is reproduced by 
summing up ladder series of diagrams, both in the Fermi 
liquid regime, and when the Fermi liquid is destroyed 
at either CDW or SDW QCP. We found, however, that 
order-by-order summation is only useful in the Fermi liq- 
uid regime, but becomes meaningless at the QCP. For the 
spin vertex, we found that Ward identity is reproduced 
near a CDW transition and near an Ising-type SDW tran- 
sition, but cannot be reproduced near a Heisenberg-type 
SDW transition. We argued that this is the consequence 
of the fact that Eliashberg theory near a Heisenberg SDW 
transition is only valid when the bosonic momentum is 
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much larger than the bosonic frequency, vfQ > ^- Ward 
identity, on the contrary, is valid in the the opposite limit 
vpq < il, where Eliashberg theory is inapplicable. 
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